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scaling for anisotropic Quantum-Critical systems 
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We propose that Kibble-Zurek scaling can be studied in optical lattices by creating geometries that 
support, Dirac, Semi-Dirac and Quadratic Band Crossings. On a Honeycomb lattice with fermions, 
as a staggered on-site potential is varied through zero, the system crosses the gapless Dirac points, 
and we show that the density of defects created scales as 1/r, where r is the inverse rate of change 
of the potential, in agreement with the Kibble-Zurek relation. We generalize the result for a passage 
through a semi-Dirac point in d dimensions, in which spectrum is linear in m parallel directions and 
quadratic in rest of the perpendicular (d — m) directions. We find that the defect density is given 
by i/T-'"''ll"-||+('*~'")''i-i where v\\,z\\ and v^,z±_ are the dynamical exponents and the correlation 
length exponents along the parallel and perpendicular directions, respectively. The scaling relations 
are also generalized to the case of non-linear quenching. 
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The Kibble-Zurek (KZ) scaling [J, i, i, i, i, @ of de- 
fect density in the final state of a quantum many-body 
system following a slow passage across a quantum crit- 
ical point, has been an exciting area of recent research. 
The KZ argument predicts that the scaling of the defect 
density is universal and is given as n -^i^vd/(nuz+i) 
where r is the inverse rate of change of a parameter, 
d is the spatial dimension and v and z are the correla- 
tion length and dynamical exponents, respectively, asso- 
ciated with the quantum critical point 0, [1] across which 
the system is swept. Following the initial predictions, 
a plethora of theoretical studies have been carried out 
0, [m, [13, El, Q [S m, [13, [II to explore the de- 
fect generation and the entropy production using differ- 



10 , quan- 
13j , along 



ent quenching schemes across critical points 
tum multicritical points (isj . gapless phases 
gapless lines [ll|, etc. The possibility of experimental 
observations in a spin-1 Bose condensate [l^ and also on 
ions trapped in optical lattices [lO, [IH has provided a 
tremendous boost to the related theoretical studies. 

Here, we propose that Kibble-Zurek scaling can be 
studied in optical lattices with fermionic atoms by cre- 
ating geometries that support Dirac, semi-Dirac and 
Quadratic Band Crossings. For example, a Honeycomb 
lattice consists of two interpenetrating triangular lat- 
tices. If the two sublattices are controlled separately, 
one can create a staggered on-sitepotential, which cre- 
ates a gap in the spectrum [l^, [l^. We will call the 
Honeycomb lattice with a staggered potential a gapped 
Graphene Hamiltonian in analogy with Graphene [22| . 
The system can be loaded with atoms when one of the 
sublattices has a much lower energy than the other. As 
the staggered potential is varied through zero, the sys- 



tem will cross through gapless Dirac point, and the cre- 
ation of defect density can be studied as a function of the 
rate of change of potential. Similarly, two interpenetrat- 
ing square-lattices can lead to Quadratic Band Crossing 



* Electronic address; ldutta@iitk. acini 

t Electronic address: |singh@physics . ucdavies . edu | 



2J] and either a deformed Honeycomb system 



2e|, 

or a 3-band system can lead to semi-Dirac points [27|, 
where the spectrum is linear along one axis and quadratic 
along another. Anisotropic quantum critical points as- 
sociated with spectra that are linear in some directions 
and quadratic in others arise in systems as diverse as 
semiconductor hctcro-structurcs 0| and Hc'^ [1^. We 
will calculate a generalized Kibble-Zurek scaling for such 
anisotropic Quantum critical systems. 

Many previous theoretical studies are on one- 
dimensional quantum spin systems some of which can 
be exactly solved via Jordan- Wigner transformation [s^l ■ 
On the other hand, the quenching dynamics of non- 
integrable spin chains have been studied using adia- 
batic perturbation theory [3| or exact diagonalization 
and time-dependent density matrix renormalization tech- 
niques [l^, [13] • Both the methods have been very suc- 
cessful in yielding exciting results for the scaling of defect 
density following a slow quantum quench in pure as well 
as random systems There are also a few results for 
the quenching dynamics of higher dimensional systems, 
e.g., in ref. Il3l| th e exact solvability of a two-dimensional 
Kitaev model [3l| was utilized. 

In this paper, we propose a kind of time-dependent 
perturbation that would be very difficult to achieve in 
solid state systems but should be possible in optical lat- 
tices. Because the underlying geometry in optical lattices 
is controlled by lasers, one should be able to change them 
in such a way that the system passes through special band 
crossings. The salient feature of our proposal on Honey- 
comb and other lattices is that the time-dependent per- 
turbation preserves crystal momentum. Thus, as long as 
the interactions between particles are weak, the system 
factorizes into independent momentum sectors. For each 
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FIG. 1: A schematic picture of a honeycomb lattice consisting 
of two interpenetrating triangular lattices shown in different 
colors. In cold-atoms, this system can be controlled by two 
sets of lasers. 



value of the momentum, there is a probability that the 
system can undergo transition from the lower to the up- 
per eigenstate. This can be analyzed using the Landau- 
Zener transition formula [s^ to provide an exact result 
for the defect creation in the quenching process. In- 
tegration over all momentum leads to the desired KZ 
scaling relations. Furthermore, quenching through the 
semi-Dirac point yields a more general form of the KZ 
scaling for spatially anisotropic quantum critical points. 
In principle, Dirac-like band crossings are of tremendous 
interest in themselves, in contexts such as topological 
insulators, [ssl [3^ But, whether one can find a perturba- 
tion that drives the system through such a crossing will 
have to be examined in individual cases. 

We focus primarily on a tight-binding model of 
fermions moving on a Honeycomb lattice. This lattice 
has two inequivalent sites per unit cell (see Fig. 1) and 
hence the system is described in the momentum space 
in terms of de-coupled 2x2 matrices [2^ . Assuming 
the system to be half-filled, the honeycomb lattice of the 
tight-binding Hamiltonian is critical at two inequivalent 
points of the Brillouin zone and there is Dirac cone of 
excitations with linear dispersion around this point. If a 
staggered on-site potential (called the mass term) which 
differentiates between two inequivalent sites is added to 
the tight binding Hamiltonian, a gap is created in the 
spectrum [2^. We assume that the magnitude of the 
mass term (denoted by fj,) is quenched linearly from — oo 
to — oo at a linear rate ^ = t/r so that the system crosses 
the gapless point /i = at a time t = Q. When /i = — oo, 
only one sublattice of the HC lattice will be occupied 
which amounts to saying that for every Fourier mode fc, 
one of the basis states is occupied. If the mass term is 
slowly changed, the state would follow adiabatically until 
the Dirac point at t = where the spectrum is gapless 
and thereafter system no longer follows the ground state 
adiabatically. We here show that the density of defect 
thus generated is in complete agreement with the KZ 
prediction. We then generalize to the case of quench- 
ing through a (d, m) semi-Dirac point [l^ at which the 
spectrum is linear in m-directions and quadratic in rest 
of the (d — m)-directions. Exploiting again the 2x2 



form of the resulting Hamiltonian and the Landau-Zcner 
transition formula we arrive at a generalized KZ scal- 
ing form l/r™''ii^il+('^-™)''i^i where v\\,z\\ and v^,z±_ 
are the dynamical exponent and the correlation length 
exponents along the parallel directions and perpendicu- 
lar directions, respectively. The same scaling relation is 
also expected to be valid for quenching through a (d, m) 
quantum LiQishitz point [35i] . 

Let us start from the general Hamiltonian of fermions 
moving on a honeycomb lattice 



H = 



k=0 



(1) 



where Hk is the reduced 2x2 matrix {Hk)ii = 
{Hk)22 = and {Hk)i2 = {Hu)\2 = -t[e*^-==" + 
2g-ik^a/2 cos(A:y-\/3a/2), where t is the hopping term and 
a is the lattice spacing. In the presence of a mass term 
that arises due to a staggered on-site potential /x [2^ , the 
reduced Hamiltonian takes the form 



fi A 
A* -n 



(2) 



where A(fc) = {Hk)i2 as defined before. The 
off-diagonal term can be simplified by expand- 
ing around the Dirac point (0, ) as A(fc) ~ 



iq.a _^ 2e-»9-°/2 C0s((:^ + Qy)^) 



1x and ky = + qy 



3 a/3 a 

y - -r qy In the limit of qx,qy 

arrive at a simpler form of Eq. [2] 



where ^ 
0, we 



Hk = 



Vf{qy + iqx) 



(3) 



where vj = 3ta/2. 

During the quenching /i = t/r, the general state vector 
at an instant t can be written as {tpkit)) = Cik{t)\lk) + 
C2k{t)\2k) where the basis vectors are |lfc) and \2k)- 
With the initial condition \Cik{t — oo)p = 1, the 
non-adiabatic transition probability due to the passage 
through the Dirac point is given as pk = \Cik{t +oo)\'^. 
Using the Hamiltonian ^ , the time evolution of the sys- 
tem is described in terms of the Schrodinger equations 



4cr.(0 
^|-^C2kit) 



-Cik{t) + Aik)C2kit) 
'-C2.k{t) + ^*ik)Cik{t). 



(4) 



The above equations describes two-time dependent levels 
zLy/{t/T)'^ + A^ approaching each other and there is an 
avoided level crossing at t = 0. The non-adiabatic tran- 
sition probability can be calculated using the Landau- 
Zener transition formula [s^ for each Fourier mode k 
given by pk = exp(— 7r| Apr). We therefore obtain 
the density of defects integrating over the modes n ~ 
j Pkd?k. The non-adiabatic transition becomes promi- 
nent only in the vicinity of the Dirac point where the 
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FIG. 2: The density of defect generated due to quenching 
through a two-dimensional Dirac point, as in Graphene, ob- 
tained through the numerical integration of the Schrodinger 
equations |4] The data clearly shows a 1/r-decay of the defect 
density which is indicated with a straight line of slope -1. 



gap closes, we therefore use the form of A valid close 
to the Dirac point given in Eq. ([6]) and extend the lim- 
its of integration from — oo to +oo. The defect density 
therefore scales as 



+ 00 



dqxdqy cxp(- 



2 / 2 



(5) 



Noting that the spectrum close to the Dirac point is 
of the form ±y^(^^ + we find that the mini- 

mum gap vanishes as /i yielding the critical exponents 
vz = 1. On the other hand, the linear dispersion at 
the Dirac point = 0) yield the dynamical exponent 
z = \. The 1/r-scaling of the defect density therefore 
is in complete agreement with the Kibble-Zurek predic- 
tion n ~ ii^yd/vz+i ^ The defect density obtained via 
numerical integration is shown in Fig. 2. 

In the optical lattices, the defect density can be mea- 
sured through the total occupation of the upper band, or 
in more detail, by the momentum distribution function, 
n(k). In our two-band Graphene system, n{k) in the 
lower-band would develop a hole near the Dirac points, 
whose size would grow with the rate of quenching. Con- 
versely, in the upper-band the particles will be concen- 
trated close to the Dirac points. The resulting momen- 
tum distribution would be highly unlike a thermal broad- 
ening. A picture of the momentum distribution function 
in the upper band is shown in an extended Brillouin zone 
in Fig. 3. 

We shall now concentrate on the quenching through 
the semi-Dirac point which provides a more general situ- 
ation of the quenching dynamics. The Hamiltonian close 
to the semi-Dirac point with a mass term can be put in 
the form [13 




FIG. 3: Momentum distribution of particles created in the up- 
per band is shown in an extended Brillouin zone. The centers 
of the six bright regions form the corners of the first Brillouin 
Zone and house the two inequivalent Dirac points. 



crosses the semi-Dirac point, the defect density scales as 

.4 1 



dqx 



(7) 



The corresponding numerical result is shown in Fig. 2. 
Generalizing to three-dimesnions [2^, where the off- 
diagonal term gets modified to A, = ro/|(7||| -I- q\/2m 



with 
given by 



the scaling of the defect density is 



dq\\dqi_ exp(-7rT(ti^gj^ 



(2m) 



T 



(8) 



The scaling relations derived in Eqs. ^ and ([8]) lead to a 
very interesting generic scaling relation of the defect den- 
sity for quenching through a (d, m) anisotropic quantum 
critical point given by the modified KZ form 



(9) 

where the critical exponents (zy||,2||) and {i'^^za_) are 
conjugate to parallel and perpendicular directions, re- 
spectively. The case with m — Q refers to the quenching 
through a band-crossing point with quadratic spectrum 
in c?-dimensions while m = d is the result for quenching 
through a Dirac point. 

We shall generalize the modified KZ scaling ^ to the 
case of non-linear variation of the mass parameter ^ = 
\t / t]"" sgn{t) . The situation is simpler as the spectrum 
IS gai aless for ^ = (i.e., t = 0) [H. Following refs. 

[la |. it is easy to show that that the probablity of 
excitation for the mode fc, pfc = |C'ifc(< — > oo)p must be 
a function of |A(fc)pT^"/"+^ where the initial condition 
If the mass term ^ is quenched as t/r and the system |Cife(t — > — cx))p = 1. Expanding around the semi-Dirac 



Hi, 



-ivf\qy \ +qi/2m 



ivf\qy\+ql/2m 



(6) 
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FIG. 4: The density of defect generated due to quenching 
through a two-dimensional semi-Dirac point obtained through 
the numerical integration. The r~^'^ ''-decay of the defect den- 
sity is indicated with a straight line of slope -3/4. 

point and extending the limits of integration form —oo 
to +00, we get 



which reduces to the form ([9]) for the linear case a = 1. 

In conclusion, we have proposed that Kibble-Zurek 
scaling can be studied in optical lattices of trapped 
fermions by creating geometries that have Dirac, semi- 
Dirac or Quadratic Band Crossings. These systems are 
typically interpenetrating lattices, which can be con- 
trolled independently. Thus, one can create a staggered 
on-site potential that leads to a gap in the spectra. 
As this staggered potential is varied through zero, one 
crosses a quantum critical point and even a slow quench- 
ing would lead to the creation of a certain density of 
defects. We have calculated the density of such defects 
for the gapped Graphene Hamiltonian and have gener- 
alized the Kibble-Zurek results to anisotropic Quantum 
Critical points that arise with semi-Dirac spectra. 

We thank Satyajit Banerjee, Debashish Chowdhury, 
Warren Pickett and Diptiman Sen for helpful comments 
and discussions. 
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